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Abstract 

States of nonlinear quantum oscillators (f-oscillators) are consid- 
ered in the Weyl-Wigner-Moyal representation and the tomographic 
probability representation, where the states are described by stan- 
dard probability distributions instead of wave functions or density 
matrices. The evolving integrals of motion for classical and quantum 
f-oscillators are found and the solution for the Liouville equation as- 
sociated with the probability distribution on the phase space for this 
oscillator is obtained along with the solution of Moyal equation for 
quantum f-oscillator, which provide the solutions for partial case of 
f-nonlinearity existing in Kerr media. Nonlinear coherent states and 
the thermodynamics of nonlinear oscillators are studied. 

keywords: coherent states, Wigner function, tomographic probability, de- 
formed oscillator. 



1 Introduction 

The small vibrations of classical and quantum systems are usually described 
by linear harmonic oscillator model. Recently [1,2] the notion of f-oscillators 
and corresponding nonlinear coherent states for the oscillator quantum do- 
main was introduced (see also [3]). The f-oscillators are nonlinear oscillators 
with a specific kind of nonlinearity for which the frequency depends on the 
oscillation energy [4]. The known q-oscillators [5,6] are very particular cases 
of the nonlinear vibrations with exponential dependence of the frequency on 
the energy of the vibrations. The properties of the f-oscillators and use of 
these oscillators for constructing models of different phenomena are inten- 
sively discussed in the literature, see, e.g., [7-18]. 
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The nonlinearity of f-oscillators provides a deformation mechanism cre- 
ating quantum group structures. For example, the f-deformed Weyl systems 
corresponding to the deformation of the Heisenberg-Weyl group were studied 
in [19] and q-nonhnearity corresponds to q-deformed Heisenberg-Weyl group. 
Recently in quantum [20] and classical mechanics [21, 22] the tomographic 
probability representation of system states was suggested. In classical sta- 
tistical mechanics, the tomographic probability distribution (tomogram or 
symplectic tomogram) is the Radon transform [23] of the probability distri- 
bution of the system on its phase space which is also the standard probability 
distribution. In quantum mechanics, the tomogram is the Radon transform 
of the Wigner function [24] of quantum states. In Moyal approach [25] , the 
quantum states are described by the Wigner function and obey the evolution 
equation which is similar to the classical Liouville equation. The Wigner 
function, which is the Weyl symbol [26] of the density operator p, can take 
negative values and, in view of this, it is not a probability distribution. It 
is called quasidistribution and contains complete information on quantum 
states. The symplectic tomogram also contains complete information on 
quantum state of a system but the tomogram is the measurable nonnegative 
probability distribution. The approach to describe the quantum states by 
the synplectic tomograms is called the tomographic probability representa- 
tion of quantum mechanics and its properties are studied in different aspects 
in [27-34]. 

The aim of this paper is to find integrals of motion and to consider the 
behaviour of f-oscillators, both classical and quantum ones, in the tomo- 
graphic probability representation. We study also quantum f-oscillators in 
the Moyal representation. We study the thermodynamics of q-oscillator in 
the case of small nonlinearity. Both the Weyl-Wigner-Moyal representation 
and the tomographic probability representation are examples of realizing the 
generic star-product scheme (see, e.g., [27,35-39]). In this context, the goal 
of the paper is also to associate the behaviour of quantum nonlinear oscil- 
lator with the star-product quantization framework. We point out that the 
electromagnetic-field vibrations in Kerr media are a particular example of 
nonlinear f-oscillations with nonlinearity function proportional to the field 
amplitude. 

The paper is organized as follows. 

In Sect. 2, the properties of f-oscillators are reviewed both in the classical 
and quantum domains. In Sect. 3, the Moyal approach and construction 
of Weyl symbols, including a deformed Wigner function, are discussed. In 
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Sect. 4, the tomography of classical and quantum states based on applying 
integral Radon transform is presented. In Sect. 5, nonlinear coherent states 
are studied. In Sect. 6, the partition function of nonlinear oscillator is 
considered. Finally the prospects and conclusions are done in Sect. 7. 



2 f-oscillators and kinetic equation 
2.1 Classical nonlinear oscillator 

The coordinate x{t) of linear harmonic classical oscillator satisfies the equa- 
tion of motion 

x + x = 0, (1) 

where for simplicity the mass and frequency of the oscillator are set to be 
equal to unity, m = u = 1. 

For complex amplitude a{t) == {l/\^){x + ix) the equation of motion 
reads 

a = —ia, (2) 

and its solutions have the form 

a{t) = aoe-'^ (3) 

with cto being any initial complex amplitude. The energy of vibrations is the 
integral of motion 

i?=Kt)r=Kr = y + y. (4) 

Thus the vibration of linear harmonic oscillator has the important property 
— the frequency of the vibration does not depend on the vibration energy. 
Another formulation of the same property implies that the phase of the vi- 
brating linear oscillator is insensitive to the absolute value of the amplitude 
of the vibration. The simplest nonlinearity of vibration just violates this 
property. Let us consider the generalization of ([3]) taking the equation of 
motion of the form 

a = —i{aa*)a. (5) 
One can see that -^(aa*) = 0, and the solutions ([3]) become 

a(t) = aoe-^l""l'*. (6) 
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Thus for cubic nonlinear equation the frequency of the vibration depends 
quadratically on the amphtude of the vibrations. The generic form of the 
equation of motion can be written as 

a = —iuj{aa*)a, (7) 

where the frequency function uj{E) depends on energy (jlj). The nonhnearity 
of vibrations is coded by the function uj{E). This means that the solution of 
Eq. © 

a{t) = aoe-'"("o"S)* (8) 

corresponds to the kind of the nonhnearity of the equation. If one describes 
the nonlinear-oscillator motion by the Hamiltonian 

H = a}af, (9) 

in which the variable «/ reads as 

«/ = af{a*a), (10) 

the dependence of the frequency on the energy is related to the dependence 
on the energy of the real function f{E). 
In fact, from (fTOj) one has 

u{E) = f{E) + Ef\E). (11) 

We call the f-oscillator the classical oscillator with the Hamiltonian ([2]), where 
a/ is given by (ITO!) . thus the oscillator evolves as 

a{t) = aoexp{-it[f (aoal) + (aoal) f'{aoal)]} . (12) 

The Liouville equation for the probability distribution f{q,p,t) on phase- 
space reads 

df{q,P,t) ^ df{q,p,t) . ^ df{q,p,t) . ^ ^ 
dt dq dp ' 

Since for our f-oscillator 

q=p[f{E)+ Ef (E)] , p=-q[f (E) + Ef (E)] , (14) 

one has the Liouville equation of the form 

dJi^^OJi^^ [/ (E) + Ef (E)]-^^, [/ (E) + Ef m = 0, 

(15) 
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where E = {q^ + p'^)/2 = {ql + pl)/2. 

The system of equations (fT^ has two time-dependent integrals of motion, 
following the solution f|T2l) of equation of motion ([7j) 



Po = p cos [u (E) t] + q sin [u (E) t] , 
go = —p sin [u (E) t] + q cos [u (E) t] . 



(16) 
(17) 



For a linear oscillator, / = 1 and the integrals of motion have the form [40] 
Po = pcost + qsint, q^ = —psint + q cost (18) 
and their dependence on complex amplitude reads 

(19) 



V2 



For f-oscillator, one has the complex integral of motion generalizing the 
above invariant 

q + ip 



V2 



exp 



t 



(20) 



In terms of the integrals of motion, the solution of Liouville equation reads 

f{q,P,t) = fo{qo{t),Po{t)), (21) 

where fo{q,p) is the probability distribution given at time t = 0. Thus, 
given the initial distribution /q (g, p) the form of the solution of the kinetic 
equation (HM is 



f{q,p,t) 



fo 



—p sm uj 



t + q cos UJ 



'q2 j^p2' 



pcosiu 



'q2 j^p2' 



+q sin u 



q2 _|_ p2 



t 



(22) 



One has the example of classical f-oscillator, the so-called classical q- 
oscillator [5,6], for which 



/ {aa* 



I sinh Xaa* 



Xaa* 



(23) 
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For small nonlinearity A ^ 1, one has the approximate value of the nonlin- 
earity function 

/ {aa*) - 1 + (aa*)^ = 1 + j^E^ E = aa* (24) 
and, in this case, the frequency of the q-oscillator is 

u{E)^l + (25) 

so it depends quadratically, for small nonlinearities, on the energy of vibra- 
tions. 



2.2 Quantum nonlinear oscillator 

Let us consider the quantum oscillator annihilation and creation operators a 
and , respectively. For Hamiltonian of the form 

H ^a)a + -M,a)\^i, (26) 

one has the solution to Heisenberg evolution equation {h — 1) 

a{t) = ae-'\ a^^{q + ip), (27) 

where q and p are the position and momentum operators, respectively. 
Let us now consider the quantum f-oscillator with the Hamiltonian 

H{n)^^ (^A\Af + AfA\) , n = a^a, (28) 

where 



Af = af (a^aj . (29) 
Now we consider the Hamiltonian of the form 

H = A\Af. (30) 

Physical phenomenon where the quantum f-oscillator Hamiltonian might 
prove to be useful is the Kerr effect (see, for instance, [41]). The description 
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of such a system has been recently considered in [42] with the Hamiltonian 
(in dimensionless units) 

Hk = X (o'^) + a) a, (31) 

where the hght frequency is equal to unity and the crystal nonlinear suscep- 
tibility is Xi with the aim of describing the evolution in terms of the Moyal 
representation. The above operator can be easily written as a function of the 
number operator h = d^a, 

Hk = h + xnih- 1) . (32) 
This suggests a description by f-oscillators introducing 

Ak = afK (n) (33) 

with 

fK (n) = ^l-X + Xn. (34) 

The function fx above introduced allows one to consider //^-nonlinear co- 
herent states and their evolution. 

The f-oscillator for both Hamiltonians ( 128|) and ( l30l) has two time-dependent 
integrals of motion of the following form: 

Q(t) = e-^^(")*ae*^(")*, Q^{t) = e-*^(")*a^e*^(")*. (35) 

The total time derivative of these operators equals zero. For time t = 0, the 
integrals of motion coincide with the annihilation and creation operators, 
respectively. 

Due to the structure of the f-oscillator Hamiltonian, the eigenstates \n, f > 
of this oscillator Hamiltonian coincide with the eigenstates of the operator 
a^a\n > related to the usual oscillator. Since the Hamiltonian is expressed in 
terms of number operator n, the result of commutation provides an explicit 
form of the integrals of motion. One can use the following formulae: 

CKf (n) = ip (n + Vj a, a^ip + = (p (n) a\ (36) 

where (p (n) is an arbitrary function. Then the Hamiltonian (128!) takes the 
form 

Hin) = ^ {hf (hf + (n + l)fih + lf). 
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The above integrals of the motion can be written as 

g(t) =aF(n,t), Q^it) = Fin,t)d\ (37) 
where the function F {n, t) for both Hamiltonians fl25]) and fl3Up reads 



F (n, t) = f (n) exp 



(38) 



[h (h) - H {h-i)}t 

The evolution of any initial density state po (ji, a^^ is given as 

p[a,a\t) =po[Q{t),Q\t)). (39) 
The above density operator satisfies the von Neumann evolution equation 

d 



dt' 



p + i H (n) ,p 



0. 



(40) 



3 Wigner functions 

The Wigner function of the evolving f-oscillator state reads 
W {q, p, t) = 2Tr [Pp {t) exp {2 (aa^ - a*aj } 



(41) 



where P is the parity operator P = (—1) and complex number a is ex- 
pressed in terms of the position and momentum as a = {q + ip)/V2. The 
parity operator commutes with the Hamiltonian. In view of this, the Wigner 
function takes the form 



W{q,p,t) =2TT[Ppoexp{2 [aQ^ (t) - a*Q{t))} 



(42) 



The Wigner function P2|) is the solution to the Moyal evolution equation [25] 
for the f-oscillator. To calculate explicitly the Wigner function, one needs to 
evaluate the trace in Eq. fl38|) . It is worth noting that in Eq. P2l) the Weyl 
system operator has the form of f-deformed shift operator which was used 
in [19]. The deformation function is given by Eq. ( 138|) . 

When the initial state is the ground oscillator state po = |0 >< 0|, one 
has the Wigner function 



Wo {q,p,t) = 2<0\ exp h (aQ^t) - a*Q{t) 



10 > . 



(43) 
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Thus the Wigner function is proportional to the expectation value of the 
deformed shift operator in the vacuum state. 

Let us now discuss the deformed creation and annihilation operators 
and Af and the corresponding deformed parity operator, more formally. 

The replacement of the usual boson operators and a by the deformed 
ones gives the deformation of the shift operator 



exp 



exp 



(ait 



a* A, 



and the deformation of parity operator 



L)"^" = exp (iTid^a^ — >■ exp (iTrA^jrAf^ 



(44) 



(45) 



Then one has the possibility of introducing deformed Wigner functions using 
f-deformation of star-product quantization scheme. Namely, one gets either 



Wi^\q,p) = 2Tr 



Ppexp{2 (ai} 



keeping the parity operator in the usual form or 



WPiq,p) 



2Tr 



exp (iirA^Af^ pexp {2 (a A\ - «*i/)} 



(46) 



(47) 



For small deformations, the Wigner functions we have introduced do not 
differ much from the standard Wigner quasidistributions. The evolution in 
time of the deformed Wigner functions can be obtained by using the integrals 
of motion (P7|) . 



4 Symplectic tomography 

Let us discuss in this section symplectic tomography [20] of classical states 
and nonlinear quantum states. 

We will focus on symplectic tomograms of f-oscillator states. 

Given a density state p, the symplectic tomogram also called tomographic 
symbol of the density operator is defined as 

w (X, fi, v) = Tr [p6 (Xi - pq-jyp)], (48) 

where X,p,u are real numbers while 1, q and p are identity, position and 
momentum operators, respectively. In terms of the Wigner function, the 
tomogram reads 

wiX,p,u)= fwiq,p)5iX-pq-up)^. (49) 
J Ztx 
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The tomogram is a probability distribution function of the position X, i.e., 

w{X,fx,iy)>0 (50) 

and 



J w (X, n, v) dX = 1. 



(51) 



The parameters fi = s cos 6* and u = sin 6' correspond to a reference frame 
in the phase space which is scaled (g sq,p s~^p) and thereafter rotated 
{sq qs cos 9 + ps~^ sin 9) . 

In classical mechanics, an analog of the tomographic symbol can be con- 
structed for an arbitrary observable A {q,p), which is a function on the phase 
space, as the Radon integral 



w_A (X, /i, z/) = J A (g, p) 5 {X — fiq — up) dq dp. 



(52) 



Let us discuss now the tomographic representation for the classical f-oscillator. 
The solution to the Liouville equation f[T^ has the following tomogram for 
the nonlinear-oscillator state: 



w (X, /i, u,t)= / /o 



-q cos uj 



' q^ + p^' 



p sm ijj 

p cos UJ 



'q2 _^p2' 



' q^ + p^' 



t 



t + qsincu 



x6 {X — fiq — up) dq dp. 



(53) 



The integral provides the solution (the tomogram) in terms of initial value 
of the tomogram wq (X, yU, u) which is given as 



Wo (X, n,u) = J fo (q, p)5{X - fxq- up) dq dp. 



(54) 



The symplectic tomogram of the f-oscillator quantum states is given in terms 
of the density operator fl39l) as 



w (X, /i, u, t) = Tr U (Qit), Q^it)) 5iX-fiq- up) 



(55) 



The quantum tomogram provides the probability distribution of the f- 
oscillator position w{x, 1, 0) = P{x) and the f-oscillator momentum w{p, 0, 1) 
Pip). 
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The ground state of the nonhnear oscillator satisfies the equation Af\0, f >— 
0, which also satisfies the equation for the usual oscillator ground state 
a|0, / >= 0. This means that the tomogram of the f-oscillator ground state 
reads 

Wo{X, //, I/) = IV [|0, / X 0, f\5{X-piq- up)] = , ^ exp ^ 



(56) 

The excited states of the nonlinear f-oscillator |n, / > have the tomogram, 
which also coincides with the tomogram of the excited state |n > of the usual 
harmonic oscillator, which reads 

1 / \ 

The constructions of Wigner functions and tomograms of of the nonlin- 
ear oscillator's quantum states uses the operators (dequantizers) which are 
ingredients of star-product quantization schemes. Thus we found the sym- 
plectic tomograms of the basis excited states for both harmonic and nonlinear 
oscillators which turn out to be the same probability distributions. 



5 f-oscillator coherent states 

In [1,2] the notion of f-oscillator nonlinear coherent state was introduced 
(see also [3]). The states are constructed as eigenstates of the operator Af — 
af (n), with fi — a^a and the equation 

Af\a, f >= a\a, f >, (58) 

where a is a complex number. For / (n) = 1, the nonlinear coherent state 
coincides with usual coherent state jo; >, i.e., 

I a >= exp 

and a^a\n >— n\n >, n = 0, 1, 2, . . . 

The nonlinear coherent state is expressed as the following series in the 
Fock number states \n > 

oo n 

\^J>=Nfj:——^\n>, /(n)! = /(0)/(l)/(2).../(n-l)/(n), 
„=o/ (^)!vn! 

(60) 



\a\ 



E 

n=0 



a 



\n > 



(59) 
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where the normahzation is given by the series 



( °° lap" \^^/^ 
Vn=o[/w!] ^7 



In the position representation, the nonhnear coherent state reads 



2 1 



(62) 



One can introduce two-mode nonUnear coherent states of the two-dimensional 
f-oscillator, for example, from operators 

Ai = oi/ {hi, hi) and A2 = 02/ (fii, ris) , (63) 

where Oi and 02 are annihilation operators for linear oscillators and hi = ajcij, 
i — 1,2. Their eigenstates are defined by the equations 

A|«i«2, / >= Oii\aia2, f > . (64) 

The nonlinearity is coded by the function of two variables / (711,712) which, 
in fact, depends on the energies of both mode vibrations. 

As an example, let us consider the nonlinearity given by the one variable 
function / (ni -|- 77,2) where the coherent states are entangled by construction. 
The explicit form of such nonlinear oscillator coherent state in the position 
representation reads 

*aia2/ [X, y)^Nf 2^ ^-j— ■ . . ■ TT < 2^^ ^11712 > ■ (65) 

The normalization constant is given by the expression 

E ' ^' ' 2 (66) 

nin2=0^l!^2! [/(ni + 7l2)!] / 



and 



g 2 2 

< xy\nin2 >= , =Hni (x) H^^ {y) . (67) 

0Fy2"i+'»2^i!7i2| 
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The functions are Hermite polynomials. For f — 1, the nonlinear coher- 
ent state becomes separable two-mode coherent state 



6 2 2 

*aia2 y) = 7^ exp 



TT 



«2 



2 + V2 (aix + a2y) -— — 



(68) 

Let us discTiss the introduced entanglement in more details. 
The pure state of two-mode system is separable if its wave function is 
factorized, namely, 

{Xl,X2) =1plixi)'tp2{x2), (69) 

where we denote the modes coordinates as Xi and X2. A wave function, which 
is the superposition of separable states, is entangled. 

The usual nondeformed two- mode coherent states are separable. In fact. 



|q!iQ!2 >= > \(^2 >— >■< a;ia;2|Q;iQ;2 >=< Xi\ai >< X2\ol2 >, 



(70) 



where 



< Xi\ai >— exp 



tti 



and 



< a;2|a2 >= exp 



\a2\ 



i2n 



E 

ni=0 



E 



ni 



(71) 



fH^, {X2) . (72) 



The nonlinear coherent states which we have introduced have the struc- 
ture 



\aia2,f>^Nf 



\n1n2 > . 



(73) 



In the case of / {rii + 722) = 1, this series is reduced to the product of two 
series, each of them depending on single variable. But this is not the case in 
presence of the nonlinear function / and therefore the nonlinearity creates 
entanglement. 
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6 Thermodynamics of f-oscillators 



For standard harmonic oscillators with Hamiltonian (!26|) . the thermodynam- 
ical properties are determined by the partition function 



Zo(/3) = Tre-^^ = ^e 

n=0 



-/3{n+l/2) 



1 



2 sinh(/5/2) 



(74) 



The other characteristics are given using the partition function, like the en- 
ergy 



E = Tt\H 



Zo if3) ) Zo {(3) 



n H — 



E 

n=0 



-/3(n+l/2) 



Zo iP) \ d/3 



(75) 



and the entropy 

S- 



-Tr 
1 



log 



Z^{(3) "Zo(/3) 



■^o(/?)S' 



-/3(n+l/2) 



/3 n + - -logZo(/3) 



which gives 



(3E + log Zo {(3), 

E-TS=-TlogZo{(3) 



(76) 
(77) 



which is the free energy. 

To consider the thermodynamics for f-oscillators, let us use the Hamil- 
tonian Hf = A^'jAf + 1/2. The partition function for such a system will be 
denoted Zf (/?). For the nonlinearity function containing a small deviation 
from the linear case, we assume the nonlinearity to be expressed as 



/(n)^(l + # (n)) 



(78) 



where g <^ 1 and (n) is some function of number of vibrations, i.e., we 
assume 



Af a[l + g^J (n)] 



(79) 
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(80) 



Then the proposed Hamiltonian of nonhnear oscillator takes the form 

Hf = A\Af + ^^d^d + ^ + gx{n), 
in which 

X{n)^2tl;{h). (81) 
Now one can calculate the partition function for this f-oscillator 

1^ 



Zf {(5) = Tr exp 



(82) 



Keeping the same accuracy, i.e., linear terms in the small parameter g, we 
get 



r / l\ 
Z^(/3) = $:exp -/5 n + - 

n=0 L ^ 



[l-Pgx{n)]. 



(83) 



Repeating the calculations as done for the linear oscillator, we get the de- 
formed partition function 



where 



Z{(3)^^ZoW) [l-(3g<x{n) >] , 
1 



,-^(n+l/2) 



(84) 



(85) 



Having the correction to the linear oscillator partition function, one can 
compute the small corrections to all the thermodynamical characteristics of 
the nonlinear f-oscillator, like free energy, entropy, etc. 
For example, for quantum q-oscillator considered in [19] 



fin) 



I sinh An 



An 



1 + T2"' 



we have 



A^ 



nl 



(86) 



(87) 



Here A is a small parameter related to g, i.e., g — A^/6. So in the present 
case, we have 



< X (^) >=< n^ > 



^ V n2e-^("+i/^) = - coth ^ ( coth ^ - 1^ 



(88) 
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Thus the correction for the partition functions for small q-oscillator nonlin- 
earity reads 

Z (/?)^ - Zo ^ coth ^ (^coth ^ - 1 j . (89) 

For large temperature /5 — > 0, the thermodynamic characteristics of the non- 
linear quantum oscillator become the thermodynamic characteristics of the 
classical q-oscillator. 



7 Conclusions 

To conclude, we point out the main results of our work. 

Wc studied nonlinear f-oscillators in both classical and quantum settings 
and found new time-dependent integrals of motion for the oscillators. For 
both settings, we obtained the solutions to the kinetic equations like Liou- 
ville classical equation for probability distribution on the phase space and 
von Neumann equation for the density operator of the f-oscillator states. Us- 
ing the deformation function determining the f-oscillator, we introduced the 
deformed Wigner functions by deforming Weyl displacement operator in the 
formalism of Moyal star-product quantization. Also we studied symplectic 
tomography of the f-oscillator states and obtained the tomographic proba- 
bility of the ground and excited states of the f-oscillators, which were shown 
to coincide with the tomograms of the excited states of the linear harmonic 
oscillator. We showed that the particular nonlinearity of vibrations we have 
considered in the case of two-mode nonlinear coherent states has created en- 
tangled states. We calculated the thermodynamic properties of the nonlinear 
quantum oscillators in the case of small nonlinearity and found explicit cor- 
rections to the partition function of the linear oscillator induced by small 
nonlinearity of q-oscillator. We showed that known Kerr nonlinearity of me- 
dia can be considered within the framework of a specific nonlinear oscillator 
formalism. 
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